## Loading required package: foreach
## 
## Attaching package: 'foreach'
## The following objects are masked from 'package:purrr':
## 
##     accumulate, when
## Loading required package: iterators
## 
## Attaching package: 'data.table'
## The following objects are masked from 'package:lubridate':
## 
##     hour, isoweek, mday, minute, month, quarter, second, wday, week,
##     yday, year
## The following objects are masked from 'package:dplyr':
## 
##     between, first, last
## The following object is masked from 'package:purrr':
## 
##     transpose
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`
## [1] 1

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.060185 -0.003642  0.000510  0.003365  0.041035 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0022339  0.0001474  -15.16   <2e-16 ***
## test_set$FCH4  0.1299140  0.0073970   17.56   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006975 on 2291 degrees of freedom
## Multiple R-squared:  0.1187, Adjusted R-squared:  0.1183 
## F-statistic: 308.5 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.170771 -0.002274  0.000127  0.001837  0.178992 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.229e-03  6.965e-05  -17.65   <2e-16 ***
## result$FCH4  5.134e-01  3.089e-03  166.22   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006624 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.7508, Adjusted R-squared:  0.7508 
## F-statistic: 2.763e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 1
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 2

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.071636 -0.003196  0.000791  0.003452  0.103163 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0026030  0.0001583  -16.44   <2e-16 ***
## test_set$FCH4  0.1006344  0.0071447   14.09   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.00754 on 2291 degrees of freedom
## Multiple R-squared:  0.07969,    Adjusted R-squared:  0.07929 
## F-statistic: 198.4 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.194724 -0.002160  0.000432  0.002121  0.085284 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.606e-03  7.047e-05  -22.79   <2e-16 ***
## result$FCH4  4.566e-01  3.125e-03  146.12   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006702 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.6996, Adjusted R-squared:  0.6995 
## F-statistic: 2.135e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 2
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 3

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.040536 -0.003163  0.000527  0.002967  0.045172 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.002112   0.000138  -15.31   <2e-16 ***
## test_set$FCH4  0.101434   0.005867   17.29   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006558 on 2291 degrees of freedom
## Multiple R-squared:  0.1154, Adjusted R-squared:  0.115 
## F-statistic: 298.9 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.161699 -0.002215  0.000296  0.001926  0.090340 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.402e-03  6.822e-05  -20.56   <2e-16 ***
## result$FCH4  4.409e-01  3.026e-03  145.71   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006488 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.6984, Adjusted R-squared:  0.6984 
## F-statistic: 2.123e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 3
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 4

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.074293 -0.003247  0.000909  0.003260  0.078962 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0028446  0.0001404  -20.27   <2e-16 ***
## test_set$FCH4  0.1094572  0.0053669   20.39   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.0067 on 2291 degrees of freedom
## Multiple R-squared:  0.1537, Adjusted R-squared:  0.1533 
## F-statistic:   416 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.168504 -0.002251  0.000502  0.002290  0.075533 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.861e-03  6.816e-05   -27.3   <2e-16 ***
## result$FCH4  4.018e-01  3.023e-03   132.9   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006482 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.6584, Adjusted R-squared:  0.6584 
## F-statistic: 1.767e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 4
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 5

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.055604 -0.003821  0.000448  0.003576  0.054082 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0025380  0.0001681  -15.09   <2e-16 ***
## test_set$FCH4  0.1029926  0.0085694   12.02   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007965 on 2291 degrees of freedom
## Multiple R-squared:  0.05931,    Adjusted R-squared:  0.0589 
## F-statistic: 144.4 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.191550 -0.002307  0.000214  0.001947  0.194081 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.309e-03  7.525e-05  -17.39   <2e-16 ***
## result$FCH4  5.135e-01  3.337e-03  153.87   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007157 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.7208, Adjusted R-squared:  0.7208 
## F-statistic: 2.368e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 5
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 6

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.059442 -0.003249  0.000606  0.003189  0.045591 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.002504   0.000144  -17.38   <2e-16 ***
## test_set$FCH4  0.107270   0.006455   16.62   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006867 on 2291 degrees of freedom
## Multiple R-squared:  0.1076, Adjusted R-squared:  0.1072 
## F-statistic: 276.2 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.111727 -0.002138  0.000367  0.002146  0.091988 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.635e-03  6.698e-05  -24.41   <2e-16 ***
## result$FCH4  4.512e-01  2.971e-03  151.88   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.00637 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.7156, Adjusted R-squared:  0.7155 
## F-statistic: 2.307e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 6
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 7

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.066273 -0.003303  0.000721  0.003176  0.037723 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0023071  0.0001374  -16.79   <2e-16 ***
## test_set$FCH4  0.1308690  0.0060818   21.52   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006534 on 2291 degrees of freedom
## Multiple R-squared:  0.1681, Adjusted R-squared:  0.1678 
## F-statistic:   463 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.176707 -0.002282  0.000344  0.002076  0.095323 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.485e-03  6.467e-05  -22.97   <2e-16 ***
## result$FCH4  4.479e-01  2.868e-03  156.18   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.00615 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.7268, Adjusted R-squared:  0.7268 
## F-statistic: 2.439e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 7
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 8

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.045952 -0.003452  0.000762  0.003332  0.061767 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0024565  0.0001495  -16.43   <2e-16 ***
## test_set$FCH4  0.1202037  0.0060581   19.84   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007116 on 2291 degrees of freedom
## Multiple R-squared:  0.1466, Adjusted R-squared:  0.1463 
## F-statistic: 393.7 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.147316 -0.002235  0.000376  0.002097  0.085438 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.544e-03  6.936e-05  -22.26   <2e-16 ***
## result$FCH4  4.350e-01  3.076e-03  141.40   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006596 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.6856, Adjusted R-squared:  0.6856 
## F-statistic: 1.999e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 8
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 9

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.043495 -0.003283  0.000498  0.003219  0.068228 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0026574  0.0001369   -19.4   <2e-16 ***
## test_set$FCH4  0.0952281  0.0052045    18.3   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006536 on 2291 degrees of freedom
## Multiple R-squared:  0.1275, Adjusted R-squared:  0.1271 
## F-statistic: 334.8 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.168388 -0.002337  0.000424  0.002329  0.080888 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -0.0017827  0.0000703  -25.36   <2e-16 ***
## result$FCH4  0.4039347  0.0031179  129.55   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006686 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.6467, Adjusted R-squared:  0.6467 
## F-statistic: 1.678e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 9
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 10

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.048676 -0.003339  0.000558  0.003213  0.108503 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0023348  0.0001481  -15.77   <2e-16 ***
## test_set$FCH4  0.1153319  0.0070928   16.26   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007024 on 2291 degrees of freedom
## Multiple R-squared:  0.1035, Adjusted R-squared:  0.1031 
## F-statistic: 264.4 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.154339 -0.002194  0.000272  0.001979  0.096263 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.392e-03  6.634e-05  -20.99   <2e-16 ***
## result$FCH4  4.698e-01  2.942e-03  159.67   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006309 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.7355, Adjusted R-squared:  0.7355 
## F-statistic: 2.55e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 10
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 11

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.058739 -0.003603  0.000674  0.003438  0.089687 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0026711  0.0001555  -17.18   <2e-16 ***
## test_set$FCH4  0.1062575  0.0083367   12.75   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007369 on 2291 degrees of freedom
## Multiple R-squared:  0.06621,    Adjusted R-squared:  0.06581 
## F-statistic: 162.5 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.160146 -0.002130  0.000347  0.001977  0.129836 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.455e-03  6.684e-05  -21.77   <2e-16 ***
## result$FCH4  4.932e-01  2.964e-03  166.37   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006357 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.7512, Adjusted R-squared:  0.7511 
## F-statistic: 2.768e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 11
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 12

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.041189 -0.003543  0.000551  0.003349  0.062674 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0024240  0.0001427  -16.99   <2e-16 ***
## test_set$FCH4  0.1085388  0.0078444   13.84   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006739 on 2291 degrees of freedom
## Multiple R-squared:  0.07712,    Adjusted R-squared:  0.07672 
## F-statistic: 191.4 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.088590 -0.002258  0.000217  0.001933  0.116428 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.279e-03  6.477e-05  -19.75   <2e-16 ***
## result$FCH4  5.009e-01  2.872e-03  174.37   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.00616 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.7683, Adjusted R-squared:  0.7683 
## F-statistic: 3.04e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 12
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 13

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.047161 -0.003142  0.000342  0.003147  0.083294 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0024662  0.0001453  -16.97   <2e-16 ***
## test_set$FCH4  0.0925520  0.0062235   14.87   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006905 on 2291 degrees of freedom
## Multiple R-squared:  0.08803,    Adjusted R-squared:  0.08764 
## F-statistic: 221.2 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.187353 -0.002254  0.000240  0.001963  0.105656 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.447e-03  7.037e-05  -20.56   <2e-16 ***
## result$FCH4  4.444e-01  3.121e-03  142.38   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006693 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.6886, Adjusted R-squared:  0.6885 
## F-statistic: 2.027e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 13
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 14

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.050861 -0.003339  0.000563  0.003200  0.084015 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0023358  0.0001497  -15.61   <2e-16 ***
## test_set$FCH4  0.1272967  0.0071459   17.81   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007116 on 2291 degrees of freedom
## Multiple R-squared:  0.1217, Adjusted R-squared:  0.1213 
## F-statistic: 317.3 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.105043 -0.002158  0.000323  0.002024  0.093241 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.494e-03  6.495e-05   -23.0   <2e-16 ***
## result$FCH4  4.616e-01  2.881e-03   160.3   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006177 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.7369, Adjusted R-squared:  0.7369 
## F-statistic: 2.568e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 14
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 15

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.065765 -0.003254  0.000696  0.003580  0.045389 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0026641  0.0001521  -17.52   <2e-16 ***
## test_set$FCH4  0.0689746  0.0065791   10.48   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007234 on 2291 degrees of freedom
## Multiple R-squared:  0.04578,    Adjusted R-squared:  0.04536 
## F-statistic: 109.9 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.200464 -0.002251  0.000304  0.002104  0.197397 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.521e-03  7.842e-05   -19.4   <2e-16 ***
## result$FCH4  4.659e-01  3.478e-03   133.9   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007458 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.6618, Adjusted R-squared:  0.6618 
## F-statistic: 1.794e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 15
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 16

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.065809 -0.003324  0.000780  0.003399  0.123005 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0027854  0.0001619  -17.21   <2e-16 ***
## test_set$FCH4  0.0989134  0.0070089   14.11   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007719 on 2291 degrees of freedom
## Multiple R-squared:  0.07998,    Adjusted R-squared:  0.07958 
## F-statistic: 199.2 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.146221 -0.002219  0.000435  0.002173  0.098295 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.684e-03  7.192e-05  -23.42   <2e-16 ***
## result$FCH4  4.527e-01  3.190e-03  141.93   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.00684 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.6872, Adjusted R-squared:  0.6872 
## F-statistic: 2.014e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 17

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.043002 -0.003407  0.000699  0.003223  0.038154 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0028823  0.0001357  -21.24   <2e-16 ***
## test_set$FCH4  0.1220999  0.0067970   17.96   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006467 on 2291 degrees of freedom
## Multiple R-squared:  0.1235, Adjusted R-squared:  0.1231 
## F-statistic: 322.7 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.173758 -0.002023  0.000428  0.002048  0.195287 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.613e-03  7.064e-05  -22.83   <2e-16 ***
## result$FCH4  5.159e-01  3.133e-03  164.69   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006718 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.7473, Adjusted R-squared:  0.7473 
## F-statistic: 2.712e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 17
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 18

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.052561 -0.003349  0.000566  0.003284  0.087575 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0026149  0.0001459  -17.93   <2e-16 ***
## test_set$FCH4  0.1001950  0.0056968   17.59   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006959 on 2291 degrees of freedom
## Multiple R-squared:  0.119,  Adjusted R-squared:  0.1186 
## F-statistic: 309.3 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.198707 -0.002229  0.000508  0.002318  0.075995 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.757e-03  7.013e-05  -25.06   <2e-16 ***
## result$FCH4  4.119e-01  3.110e-03  132.44   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006669 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.6567, Adjusted R-squared:  0.6567 
## F-statistic: 1.754e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 18
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 19

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.053435 -0.002855  0.000555  0.002953  0.094091 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0022447  0.0001495  -15.01   <2e-16 ***
## test_set$FCH4  0.1145451  0.0057200   20.02   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007128 on 2291 degrees of freedom
## Multiple R-squared:  0.149,  Adjusted R-squared:  0.1486 
## F-statistic:   401 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.175257 -0.002265  0.000385  0.002136  0.083069 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.582e-03  7.014e-05  -22.56   <2e-16 ***
## result$FCH4  4.157e-01  3.111e-03  133.64   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.00667 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.6608, Adjusted R-squared:  0.6607 
## F-statistic: 1.786e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 19
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 20

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.042982 -0.003262  0.000772  0.003380  0.107995 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0028766  0.0001485  -19.38   <2e-16 ***
## test_set$FCH4  0.1067937  0.0064476   16.56   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.00708 on 2291 degrees of freedom
## Multiple R-squared:  0.1069, Adjusted R-squared:  0.1066 
## F-statistic: 274.3 on 1 and 2291 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.190505 -0.002153  0.000478  0.002180  0.078782 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -0.0017428  0.0000693  -25.15   <2e-16 ***
## result$FCH4  0.4487030  0.0030735  145.99   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006591 on 9169 degrees of freedom
##   (18521 observations deleted due to missingness)
## Multiple R-squared:  0.6992, Adjusted R-squared:  0.6992 
## F-statistic: 2.131e+04 on 1 and 9169 DF,  p-value: < 2.2e-16
## 
## [1] 20
## Warning: `as.tibble()` was deprecated in tibble 2.0.0.
## ℹ Please use `as_tibble()` instead.
## ℹ The signature and semantics have changed, see `?as_tibble`.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
##    StartDate             EndDate AvailNum    EffNum ObsNum AnnNEE_gCO2
## 1 2021-09-04 2021-12-31 23:30:00     5712  535.3052   3034    464.6739
## 2 2022-01-01 2022-12-31 23:30:00    17520 1206.8784   8440   -502.5467
## 3 2023-01-01 2023-04-03 11:30:00     4440 1437.5766   2133    583.4427
##   UncNEE_Rand_gCO2 UncNEE_uStar_gCO2 UncNEE_Comb_gCO2 AnnNEE_gC UncNEE_Rand_gC
## 1        42.636103          7.981962         43.37683  126.8166       11.63604
## 2        84.754166         33.961872         91.30541 -137.1526       23.13070
## 3         8.556503          8.146882         11.81463  159.2304        2.33520
##   UncNEE_uStar_gC UncNEE_Comb_gC
## 1        2.178399      11.838197
## 2        9.268713      24.918638
## 3        2.223408       3.224392
##    StartDate             EndDate AvailNum    EffNum ObsNum AnnFCH4_gCH4
## 1 2021-09-04 2021-12-31 23:30:00     5712  901.5999    965   -0.1012138
## 2 2022-01-01 2022-12-31 23:30:00    17520 4534.3640   6122   -0.2989308
## 3 2023-01-01 2023-04-03 11:30:00     4440 1201.5861   1480   -0.1774757
##   UncFCH4_Rand_gCH4 UncFCH4_uStar_gCH4 UncFCH4_Comb_gCH4  AnnFCH4_gC
## 1        0.14065359         0.13619315        0.19578561 -0.07579047
## 2        0.14441551         0.07587745        0.16313561 -0.22384401
## 3        0.07083016         0.02912431        0.07658419 -0.13289658
##   UncFCH4_Rand_gC UncFCH4_uStar_gC UncFCH4_Comb_gC
## 1      0.10532358       0.10198353      0.14660729
## 2      0.10814057       0.05681821      0.12215847
## 3      0.05303872       0.02180874      0.05734742